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Abstract. We apply the symmetry reduction method of Roberts to numerically analyze 
the linear stability of a one-parameter family of symmetric periodic orbits with regulariz- 
able simultaneous binary collisions in the planar pairwise symmetric four-body problem 
with a mass m G (0, 1] as the parameter This reduces the linear stability analysis to the 
computation of two eigenvalues of a 3 x 3 matrix for each m e (0, 1] obtained from numer- 
ical integration of the linearized regularized equations along only the first one-eighth of 
each regularized periodic orbit. The results are that the family of symmetric periodic orbits 
with regularizable simultaneous binary collisions changes its linear stability type several 
times as m varies over (0, 1], with linear instability for m close or equal to 0.01, and linear 
stability for m close or equal to 1 . 



1. Introduction 

In Hamiltonian systems like the Newtonian A^-body problem, linear stability of a peri- 
odic orbit is necessary but insufficient for its nonlinear stability [6 1. When the periodic orbit 
is not a relative equilibrium, the characteristic multipliers are typically found by computing 
its monodromy matrix, i.e., by numerically integrating the linearized equations along the 
periodic orbit over a full period (in which the periodic orbit and its period are typically 
computed numerically as well). For a symmetric periodic orbit, Roberts ifTOl developed 
a symmetry reduction method by which the nontrivial characteristic multipliers are com- 
puted by numerical integration of the linearized equations along the periodic orbit over a 
fraction of the full period. He applied this symmetry reduction method to show that nu- 
merically the Montgomery-Chenciner figure-eight periodic orbit with equal masses f3\ is 
linearly stable 1 10 1; the numerical integration of the linearized equations along the periodic 
orbit only needed to go over one-twelfth of the full period. 

We apply Roberts' symmetry reduction method to a one-parameter family of symmet- 
ric singular periodic orbits in the planar pairwise symmetric four-body problem (PPS4BP) 
where the parameter is a mass m E (0, 1] and the singularities are regularizable simulta- 
neous binary collisions (SBCs). We recall in Section |2] the notation we used in 121 for 
the PPS4BP. (The PPS4BP is the Caledonian symmetric four-body problem | lH without 
its coUinear restrictions on the initial conditions.) To compute the nontrivial character- 
istic multipliers of these periodic orbits we numerically integrated the linearized regular- 
ized equations along each regularized periodic orbit over only one-eighth of its period. 
This shows that numerically these symmetric singular periodic orbits experience several 
changes in their linear stability type (linearly stable, spectrally stable, or linearly unstable) 
as m is varied over (0, 1]. 
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This is a marked improvement over our previous numerical investigations of the linear 
stability of these symmetric singular periodic orbits. We numerically computed IS] the 
monodromy matrix and its eigenvalues for each regularized symmetric periodic orbit start- 
ing at m = 1.00 and decreasing by 0.01 until m = 0.01. This seemed to indicate that the 
periodic orbits were linearly stable for m in the interval [0.54, 1.00] and linearly unstable 
for m in the interval [0.01,0.53]. This agreed with the stability and instability suggested 
by our long-term numerical integrations of the regularized equations starting at a numer- 
ically computed approximation of each periodic orbit's initial conditions. However, our 
numerical estimates of the monodromy matrices failed to accurately account for the trivial 
characteristic multiplier 1 of algebraic multiplicity 4: instead of getting 1 as an eigenvalue 
for each monodromy matrix, we were getting two pairs of eigenvalues, one pair of positive 
eigenvalues with one larger than 1 and the other smaller than 1 , and one pair of complex 
conjugates close to 1. As m passed below 0.61, the real eigenvalues began to move away 
from 1, so much so, that below m ~ 0.21, we had a real positive eigenvalue of the mon- 
odromy matrix whose value exceeded the limits of MATLAB. This calls into question the 
conclusions of our first attempt at determining for what values of m the symmetric singular 
periodic orbits in the PPS4BP were linearly stable and linearly unstable. 

We thus proceed to use Roberts' symmetry reduction method because it factors out, in 
an analytic manner, two of the trivial characteristic multipliers, leaving the numerical com- 
putations to estimate the two pairs of nontrivial characteristic multipliers and one pair of 
trivial characteristic multipliers. The details of these computations are given in Section]?] 
Two surprises here are the intervals [0.21,0.22] and [0.23,0.26] where we have linear sta- 
biUty. Our long-term numerical integrations of the regularized equations for these periodic 
orbits (starting at our numerical approximations of their initial conditions and over 100932 
periods) suggested instability for m in these two intervals. We also refined the numerical 
computation for m between 0.53 and 0.54 by increments of 0.001 to get a better estimate 
of that value of m where the linear stability type changes. This showed that we have linear 
stability for m — 0.539 and linear instability for m in [0.531,0.538]. 

Such changes in the linear stability type of mass-parameterized families of symmetric 
periodic orbits with regularizable collisions have been found in other A^-body problems. 
The Schubart orbit in the collinear three-body problem lfT6]| . ||4l, Q, li20J . I.17J has the in- 
ner body alternating between binary collisions with the two outer bodies. These are linearly 
stable for certain choices of the three masses f5\. Linearly stable non-Schubart orbits have 
also been found in the collinear three-body problem for certain choices of the masses 1 1 1 .| , 
lfT2ll . lfT3l . The Schubart-like orbit in the collinear symmetric four-body problem ifTSl . |fT9l . 
(TTl, f9l, fTTl, alternates between a binary collision of the two inner bodies and a SBC of 
the two outer pairs of bodies. If the masses in the collinear symmetric four-body problem 
are, from left to right, 1, m, m, and 1, then linear stability occurs when < m < 2.83 and 
m > 35.4 with linear instability for 2.83 <m < 35.4 by numerical computation of their lin- 
ear stability indices 1 19 J (a method which requires numerical integration of the regularized 
equations over a full period) and corroborated by Roberts' symmetry reduction method 
Illl. The symmetric singular periodic orbit in the fully symmetric planar four-body equal 
mass problem |8|, ifTTl (in which the position of one of the bodies determines the positions 
of the remaining three bodies) alternates between distinct SBCs and has been shown to be 
linearly stable, with respect to symmetrically constrained linear perturbations, by Roberts' 
symmetry reduction method |1 1. The linearly stable symmetric singular periodic SBC or- 
bit with m = 1 in the PPS4BP is the analytic extension |2| of the linearly stable symmetric 
singular periodic orbit in the fully symmetric planar four-body equal mass problem 1 1 1. 
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2. THEPPS4BP 

We recall from |2J the relevant notations for the PPS4BP, its regularized Hamiltonian, 
and properties of the regularized one-parameter family of symmetric SBC period orbits. In 
the PPS4BP the positions of the four planar bodies are 

{xi,X2), (X3,X4), (-Xi,-JC2), (-JC3,-X4), 

where the corresponding masses are l,m, 1, m with <m < 1. With f as the time variable 
and ■ = d/dt, the momenta for the four bodies are 

(c0i,C02) =2(ii,i2), (033,034) =2m(i3,i4), -(031,0)2), -(0)3,034). 

The Hamiltonian for the PPS4BP is 

J[032 + 03^]+^[o3? + O3i] 
1 2m 

2^x\+xj yyix3-xi)2 + {x4-X2)^ 

2m m^ 

V/(X1+X3)2 + (X2+X4)2 2^xl+xl 

The angular momentum for the PPS4BP is 

A — Xi032 — X2 03i +X3O34 — X4O33. 

A regularizable simultaneous binary collision occurs when xj, — xi ^ and X4 = X2 ^ 
(in the first and third quadrants), and also when xj — —x\ ^ and X4 = —X2 7^ (in the 
second and fourth quadrants). Initial conditions for the symmetric SBC periodic orbits in 
the PPS4BP when m = 1 are given in |2 1, and when m = 0.539 they are 

XI =2.11421, X2 = 0, X3 =0, X4 = 1.01146, 
03i =0, 0)2=0.18151, 033 =0.70392, 034 = 0. 

2.1. The Regularized Hamiltonian. We define new variables ui, U2, uj, U4, vi, V2, v^,, and 
V4 related to the variables xi,X2, X3, X4, 03i, 032, 033, and 034 by the canonical transformation 

Xl = (1 /2)(mi — M2 + M3 — M4) 

X2 ~ U\U2 + M3M4, 

X3 = (1 /2)(m3 — M4 — + M2), 
X4 = M3M4 — M1M2, 

_ V'lMl -V2M2 + V1M2-I-V1M2+V2M1 

2(mj+M2) ' 



y3M3 — V4M4 + V'3M4 + y4M3 



2(Mf + 


ul) 


7 


— V'lMl + V2M2 


+ VIU2 - 


hV2Ul 




huj) 




— V3M3 +V4U4 


+ V3U4 - 


\-v4u3 


2(1*3 " 


f «!) 
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In extended phase space, the variables are mi, M2, M3, M4, E, vi, V2, ^'3, V4, and f, where £ is 
the energy. If we set 



Ml = V'lMl - V2M2, 

M3 = V3M3 — V4M4, 



Ms = u\ — U2- 
M-i = u\ — u\- 



2 2 
- M3 M4, 



^4, 



M2 = V1M2 + V2M1, 
M4 = V3M4 + V4M3, 
Mg = 2miM2 +2M3M4, 
Mg = 2u\U2 — 2M3M4, 



then the regularized Hamiltonian for the PPS4BP in extended phase space is 



dt 
ds 



(H-E) 



1 

16 



1 + - (vf + vi)iui + ui) + {vi + vi){ui + 4) 



(M3M1 +M4M2) 



2 2 
- M2 + M3 



m {u\ + M2)("3 + "4) 



M^- 



-M2 



M4) 

M2)(«3 + «4)) 



where 



dt 

ds 



{u\ + u\){u\ + u\) 



is the regularizing change of time for this Levi-Civita regularization. The angular momen- 
tum in the new variables is 



A = - [ — V\U2 + V2U\ — V3M4 + V'4M3] . 



Let ' = d/ds. 



J = 



/ 

-/ 



for / the 4x4 identity matrix, and V be the gradient with respect to the variables 

Z= (mi,M2,M3,M4,Vi,V2,V3,V4)- 

The regularized Hamiltonian system of equations with Hamiltonian Y is 
(1) z' = /Vf(z). 

The energy E is conserved because 

dt 

E'= — =0. 
dt 

2.2. The Symmetric Periodic SBC Orbits in the Regularized PPS4BP. For m = I, 

we have analytically proven |2| the existence and symmetries of a symmetric periodic 
SBC orbit y(.s; 1), with period T = 271, £ « -2.818584789, and A = for the regularized 
PPS4BP on the level set F = 0. The initial conditions of y(i; 1 ) at i = satisfy 

M3(0;l) =Mi(0;l), m4(0;1) = -m2(0;1), 
y3(0;l) = -vi(0;l), y4(0; 1) = y2(0; 1). 
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The symmetries of "/(.s; 1) are 5f y(i; 1) = y{s + iz/2; 1) and Scjis; 1) = y(7t/2 — i; 1) where 

-GOO " 
G 
G ' 

-G 

for 



F = 



"-1 0' 




'I 


, G = 


1 




1 



Using a scaling of periodic orbits in the PPS4BP, we Q numerically continued the sym- 
metric SBC periodic orbit y(i; 1) to symmetric periodic SBC orbits 7(5; m) with A = 
for < m < 1 at 0.01 decrements with fixed period T = 2k and varying energies E{m) 
using trigonometric polynomial approximations that ensured the symmetries Sf'Y{s; m) — 
y(i + %/2;m) and Sajis; m) = 7(71 /2 — s;m). For all <m < 1 , the components of y(0; m) 
satisfy 

(2) M3(0;m) = Ml (0;m), U4{0;m) — ~U2{0;m), 

(3) V3(0;m) — — vi(0;m), V4(0;m) = V2(0;m). 

For all < m < 1, regularized SBCs occur at i = k/4,3k/4,5k/4,7k/4, where at the 
first and third times we have V3 + V4 = while at the second and fourth times we have 
v\ + = 0. The regularized symmetric periodic orbit j{s;m), in going from ^ = to 
s = 2k, corresponds in the original Hamiltonian system in the physical plane to two full 
periods of oscillation of a symmetric singular periodic orbit, whose only singularities are 
regularizable SBCs. 

Each regularized symmetric periodic orbit "/{s; m) has the trivial characteristic multiplier 
1 of algebraic multiplicity at least 4. This is because the regularized Hamiltonian Y and 
the angular momentum A are first integrals for the regularized Hamiltonian system ([T]i, and 
because of the time translation along the periodics orbit and SO(2) rotations of the periodic 
orbits (see I16J). 

3. Linear Stability of Periodic SBC Orbits 

We apply Roberts' symmetry reduction method (W\ to the one-parameter family of 
periodic orbits 'Y{s;m), Q<m< 1, of fixed period 2k, in the regularized Hamiltonian system 
M. Let V^f denote the symmetric matrix of second-order partials of f with respect to the 
components of z- It is easily shown, that if Y{t) is the fundamental matrix solution of the 
linearized equations along y(i;m), 

for an invertible Yq, then the eigenvalues of Yq^Y{2k) are indeed the characteristic multi- 
pliers of y( 

3.L Stability Reductions using Symmetries. We use the symmetries of y(i;m) to show 
that Yi^^Y{2k) can be factored in part by terms of the form Y(k/4), that is, one-eighth of 
the period of y(i;m). Thus the symmetries of y(i;m) will reduce the analysis of its linear 
stability type to the numerical computation of F (7i/4). 

Lemma 3.1. For each <m < 1, there exists a matrix W such that Yq^Y{2k) — where 
W = AD for involutions A and D with A = Y^^S^-SgYq and D = B'^ScBfor B = Y{k/4). 






F 








F 




















F 








-F 






Sg = 
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Proof. Each y{s;m) satisfies SFy{s;m) — "/{s + K/2;m). Then (by ifTOl . see also |[T]), we 
have that 

Y{kn/2)^S'pYo{Y,-'s^pY{n/2)f 
holds for all e N. Since Sp = /, taking k = 4- gives 

(4) Y{2n)=YoiYf;'s]^Y{n/2))\ 

Furthermore, each 'Y{s;m) satisfies Scjis^m) — y(7i/2 — s;m). Then (by ifTOI . see also HI), 
for 

B = Y{k/A) 

we have that 

(5) y(7t/2) = SgYoB-^SIb = SgYqB-^SgB, 

where we have used Sq = Sg- Combining equations (|4ji and (jsjl gives the factorization 

Y{2%)=YoiYf^^SlSGYoB-'SGB)\ 

By setting 

Q = S^Sg and W = Yq-^QYqB-^SgB, 

we obtain 

Fo-iF(27t) = (Y^-'QYoB-'SgB)^ = w\ 

where 



A = Y^^ ' gFo and D = ' SgB 



are both involutions, i.e., A = D 



□ 



3.2. A Choice of Yq. The matrix Q — S^Sg that appears in A is orthogonal since Sf and 
Sg are both orthogonal. Furthermore, Q is symmetric and its eigenvalues are ±1, each of 
multiplicity 4. An orthogonal basis for the eigenspace ker(2 — /) is 


















1 












-1 






















1 












1 












-1 












1 






















1 












1 


















and an orthogonal basis for the eigenspace ker(2 + /) is 


















-1 







1 



























1 







1 












1 



























-1 







1 



























1 








We look for an appropriate choice of io such that 

(6) A = Yo-'QYo-- 
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Lemma 3.2. There exists an orthogonal and symplectic Fq such that Equation Q holds. 

Proof. Since the components of 7(5; m) satisfy the Equations (j2| and (|3j, then using the 
Hamiltonian system (hj) on the level set f = 0, the components of y'(0;m) satisfy 

Uj^{Q;m) = — Mj(0;ot), M4(0;m) = M2(0;m), \'j^{Q\m) = v''i(0;m), \\{Q);m) = — v'2(0;m). 

It is easily recognized that the vector y'{Q;m) belongs to ker(2 + Now set 



and define Yq by 



(7) 



= m2(0;ot), 


c — 


v''i(0; 


m), d 


= v4(0; 


m), 


e = 


c d 


a 


b 


a 


-b 


— c 


d 


d ~c 


b 


—a 


b 


a 


~d 


~c 


c d 


a 


b 


—a 


b 


c 


~d 


-d c 


-b 


a 


b 


a 


-d 


—c 


—a b 


c 


-d 


c 


d 


a 


b 


—b —a 


d 


c 


d 


— c 


b 


—a 


a —b 


— c 


d 


c 


d 


a 


b 


—b —a 


d 


c 


~d 


c 




a 



Let col,(}o) denote the /* column of Yq. Notice that col5(Fo) = y'iO;m)/\\j'{0;m)\\. The 
last four columns of Yq form an orthonormal basis for ker(2 + /), while the first four 
columns of Yq form an orthonormal basis for ker(2 — /). Since Q is symmetric, its two 
eigenspaces are orthogonal, and so Yq is orthogonal. Note that 7col4+,(}o) = col,(Fo) for 
i = 1,2,3,4; in other words, multiplication by J maps ker(2 — /) bijectively to ker(2 + /)- 
For Pi the lower right 4x4 submatrix of Yq and P2 the upper right 4x4 submatrix of io, 
we have 



where PT" Pi +P2^P2 



Yo = 



I and Pf P2 





Pi 




Pi 


)- 


■ Pi Pi 


Pl_ 




Pl_ 




-Pi Pl_ 



0. These implies that Yq is symplectic. 



□ 



3.3. The Existence of K. By Lemma|3.1|we have Yq^Y{2%)= where W = AD with 

there exists an orthogonal 



3.2 



A = Yq^QYq cLud D = B-^SgB for B ^Y {%/4). By Lemma : 
and symplectic Yq such that Equation (|6]) holds. Choose Yq as given in Equation (|7]i. The 
matrix W = AD is then symplectic, i.e., W^JW = J, because A is symplectic with mul- 
tiplier — 1, A^JA = — J, and Sq is symplectic with multiplier —1, S^JSq = —J, and B is 
symplectic. 

Lemma 3.3. With the given choice o/Yq, there exists a matrix K uniquely determined by 
B = y(7t/4) such that 







Proof. Since W = AD where A and D are involutions, it follows that 

W"' =DA. 

By the choice of Yq, the form of the matrix A is given in Equation ([6]l. If we partition the 
symplectic matrix B into the four 4x4 submatrices. 



(8) 



B 



A3 A4 
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then the form of the inverse of B is 



^4 



^2 



Set 



Then we have that 



H 



A3 A[ 



G 
G 



D = B^SgB = 



K' Li 

-L2 -K 



where K^AIhA2 +A[i/A4, Li = A[//A2 +AIhA^, and L2 = A^//Ai +A[//A3. It fol- 
lows that K is uniquely determined by B, that 



(9) 

and that 

Thus 
(10) 



W =AD 



7 


0' 






Li' 






Li 





-/ 




.^2 


-K 




.^2 


K 



DA 









I ■ 








M 


-K 




-/ 






K 



K'^ 
K 



for a K uniquely determined by B = Y{ti/A) as was desired. 



□ 



It has been shown ifTOl that the symplectic matrix W is spectrally stable, i.e., all of 
its eigenvalues have modulus 1, if and only if all of the eigenvalues of K are real and 
have absolute value smaller than or equal to 1. The particular relationship between the 



eigenvalues of W and K given tacitly in Lemma 3.3 is as follows. The map f :C^C given 
by f(k) = (1 /2)(k+l /A.) takes an eigenvalue of W to an eigenvalue of (1/2)(W 
Note that the map / satisfies /(A.) = f{l/X). For an eigenvalue A- of W, the eigenvalue 
f(k) of (1/2)(W + W^^) is an eigenvalue of K. If A is an eigenvalue of the symplectic 
matrix W, then 1 /A, A, and 1 /A, are also eigenvalues of W. When A has modulus one, then 
A = 1/A, and 1/A = A,, and so /(A) = /(A-) which is a real number with absolute value 
smaller than or equal to 1. Thus a complex conjugate pair of eigenvalues of W of modulus 
one corresponds to a real eigenvalue of K with absolute value smaller than or equal to 1. 
When A is real, it is nonzero because W is symplectic, and /(A) = /(I /A) which is a real 
number with absolute value greater than 1 . Thus a reciprocal pair A and 1 /A of real nonzero 
eigenvalues of W corresponds to a real eigenvalue of K with absolute value greater than 1. 
When A is not real and has a modulus other than 1, then /(A) —f{l /A) and /(A,) = /(1/A-), 
with /(A) and /(A-) as complex conjugate eigenvalues of K with nonzero imaginary part. 
Thus, the four eigenvalues A, 1 /A, A, and l/XofW correspond to a complex conjugate pair 
of eigenvalues of K with nonzero imaginary part. 



3.4. The Form of K. We will show that one of the eigenvalues of is 1, and the remaining 
three eigenvalues of K are determined by the lower right 3x3 submatrix of K. Let c, denote 
the column of B = Y(n/4). 
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Lemma 3.4. With the given choice o/Yq, the matrix K uniquely determined byB = Y{%/A) 
is 

1 * * * 

c\SgJc(, C2SgJci c^SgJc^ 

c^SgJce c^SgJc-i c^ScJcfi 

c^ScJce cISgJci c\SgJc%_ 

Proof. We begin by showing that 1 is an eigenvalue of W by identifying a corresponding 
eigenvector Since Y{ti/2) — SgYoB^^SgB (Equation|5jl and Q = SpSc, it follows that 

W =Y^:^^QYoB-^SgB 

= Yq-^sISgYoB-'SgB 



Set 



■ Y,-'s]:Y{n/2). 



The orthogonality of Yq and col5(Fo) = Y {0',m) /\\y' {0;m)\\ imply that 

v = }'oV(0;m) = ||y'(0;m)||e5, 

where es = [0,0,0,0,1,0,0,0]^. Since Y{s) is a fundamental matrix, then y'{s;m) 
y(s)Fo"'Y'(0;m). Hence, 



Wv = Y^^^SIy {Ti/2)v 

,-lcT 



= Y^;'S'pY{iz/2)Y^^i{Q-m) 
= Y^'sli{%/2-m). 
Since SFj{s;m) = j{s + %/2;m) and Sp^ ~ S]^, we have that 

Y{s;m)=Sp^Y{s + n/2;m)=S'^Y{s + n/2; 



Setting i = in this gives 
From this it follows that 



Y{0;m)^Si-Y{'^/2;m 



Wv^Yf^'S'pY(.T^/2;'n) 

= yo"V(0;m) 



Thus 1 is an eigenvalue ofW and v — ||Y'(0;m)||e5 is a corresponding eigenvector 

Next, we show that the first column of K is [1,0,0,0]^. Since Wv = v, then Wes — ej. 
From the form of W given in Equation (j9]l, it follows that 



This implies that 



es = Wes 



K 



ii [1,0, 0,0" 
.^[1,0,0,0] 



"1" 




'I 



























from which it follows that the first column of is [1,0,0,0] 
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Finally we show that the lower right 3x3 submatrix of K has the prescribed entries. 
Since Yq is symplectic, the matrix B = Y{%/A) is symplectic. Hence B satisfies J = B^JB, 
and so 



B 



-JB^J. 



For W = AD with D = B SqB where Sg satisfies SqJ — —JSg we then obtain 

W ^AB-^SgB 
^A{-JB^J)SgB 
= -AJB^JSgB 
= -AJB"^ {-SgJ)B 
= AJB^SgJB. 

Writing B in the block partition form given in Equation ([8]), it follows that 



(11) 



AJB^ = 



"0 


/■ 


B^ = 





/■ 




A[ 


At 








/ 





/ 







Aj 


Al_ 




.^1 


Al 



Let co1,(5g/B) denote the /* column of SgJB. Then co1,(5g/B) — SgJci where c,- is the 
column of B = F(7i/4). This and Equation (111 imply that the (ij) entry of W is 



then cJSgJc]. But Equation (joj implies that the (6,6) entry of is the (2,2) entry of K. 
Continuing in this manner we find the remaining entries of the lower right 3x3 submatrix 
of K to be given as prescribed. □ 

3.5. A Stability Theorem. The characteristic multipliers of y(i;m) are the eigenvalues of 
which are the fourth powers of the eigenvalues of W. As was shown in the proof of 



Lemma 3 .4 an eigenvalue of /iT is L Because of Equation ( 1 1, an eigenvalue of W is 1 with 
algebraic multiplicity (at least ) 2. This accounts for two of the four known eigenvalues of 
1 for W*. Our numerical calculations show that —1 is an eigenvalue of K and hence of W 
for all <m < 1. This accounts for the remaining two known eigenvalues of 1 for W^. 

When W is spectrally stable, the eigenvalues of K are the real parts of the eigenvalues of 
W. If is an eigenvalue of K, then ±/ are eigenvalues of W and so the algebraic multiplicity 
of 1 as an eigenvalue of W'^ is at least 6. If 1 /a/2 is an eigenvalue of K, then 1 /v^ ± '/ \/2 
are eigenvalues of W, and if —l/Vl is and eigenvalue of K, then — 1/-\/2± i/v^ are 
eigenvalues of W; both these imply that —1 is a repeated eigenvalue of W^. So when the 
remaining two eigenvalues A-i and A-2 of K are real, distinct, have absolute value strictly 
smaller than one, and none of them are equal to or ± 1 / \/2, then the symmetric periodic 
SBC orbit is linearly stable, i.e., W, and hence W'^, is spectrally stable as well as semisimple 
when restricted to the four dimensional W-invariant subspace of determined by the 
two distinct modulus one complex conjugate pairs of eigenvalues of W. On the other 
hand, if one of A-i or A,2 is real with absolute value bigger than 1, or is complex with a 
nonzero imaginary part, then the symmetric periodic SBC orbit is not spectrally stable, but 
is linearly unstable. The proof of the following result about the linear stability type for 
the symmetric periodic SBC orbits in the PPS4BP follows from all of the Lemmas and 
subsequent comments presented in this Section. 

Theorem 3.5. The symmetric periodic SBC orbit j(s;m) of period T = 2% and energy 
E{m) is spectrally stable in the PPS4BP if and only ifXi and Xi ore real and have absolute 
value smaller or equal to I. IfXi and A-2 are real, distinct, have absolute value strictly 
smaller than 1, and none of them are equal to or ±1 / \/2, then 'Y{s;m) is linearly stable 
in the PPS4BP. 
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Figure 1. The eigenvalues A-i, A,2, A,3, when real, of the 3x3 lower 
right submatrix of K over < m < 1 . 



4. Numerical Results 

We computed Y{%/A) using our trigonometric polynomial approximations of 'Y{s;m) 
for each m starting at m — I and decreasing by 0.01 until we reached m = 0.01, and the 
Runge-Kutta order 4 algorithm coded in MATLAB, with a fixed time step of 

n/4 n 
50000 ^ 200000' 

From the needed columns of Y{k/4-), we computed the entries of the lower right 3x3 
submatrix of K as given in Lemma [34] and then computed the eigenvalues A-i, A,2, and X3 
of this 3x3 matrix. We have plotted these three eigenvalues, when real, as functions of m 
in Figure[T| One of these eigenvalues is real and stays close to —1 for all m e (0, 1] except 
at m — 0.20; label this eigenvalue A,3. 

The remaining two eigenvalues A-i and X2 of K that determine the linear stability type 
of j{s;m) are for m = 0.01 near 1 and not shown, respectively, in Figure [T| The values 
of A,i and X2 at m = 0.01 ai-e 0.9743145796, and -50.70044516 respectively. As m in- 
creases from 0.01, the value of A-i decreases, crossing l/\/2 for some m in (0.09,0.10), 
and crossing for some m in (0.26,0.27), while A,2 increases to the value —1.146019443 at 
m = 0.19, momentarily disappearing at m — 0.20, reappearing at m = 0.21 with a value of 
—0.8641436215, continuing to increase, crossing — 1/\/2 for a value of m in (0.22,0.23), 
until at some value of m in (0.26,0.27), we have A-i = A,2 < 0. For m in [0.27,0.53], the 
eigenvalues A-i and X2 form a complex conjugate pair with nonzero imaginary part, and 
thus disappear in Figure [T| For some value of m in (0.53,0.54), we have Xi and X2 reap- 
pearing in Figure [T] with A-i = A-2 > 0. As m increases from there, A-i increases and X2 
decreases, with Xi crossing for a value of m in (0.54, 0.55), and with the values of Xi and 
X2 at m — 1 being respectively, 



(12) 



0.6941364299, -0.6802222699, 
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where the first of these is slightly smaller than l/\/2, and the latter is slighter larger than 
— 1 / \/2- These changes in the values of A-i and A,2 account for the changes in the linear 
stability type of y{s;m) as m varies over (0, 1] 



From the numerical results and Theorem 3.5 we conclude that the periodic orbit y(i; m) 



is linearly stable when m is in [0.21,0.22], or m is in [0.23,0.25], m ~ 0.54, or m is in 
[0.55, 1]. We have linear instability when m is in [0.01,0.19] or in [0.27,0.53]. We have at 
least spectral stability when m — 0.20 where A,3 disappears momentarily along with A,2 to 
form the complex conjugate pair with nonzero imaginary part, 

-0.9972588720 ±0.008650400165/. 

This appears numerically to be a repeated eigenvalue of — 1 for K. We also have at least 
spectral stability for a value of m in (0.22,0.23), and for a value of m in (0.54,0.55). 

We have confirmed that numerically the equal mass symmetric periodic SBC orbit 
y(i; 1) is linearly stable in the PPS4BP. From the eigenvalues of K, which are 1, —1, and 



those listed in ( 12 1, the characteristic multipliers of y{s; 1) are 1 with algebraic multiplicity 
4, and the two distinct complex conjugate pairs of modulus one, 

-0.9888710746 ±0.1487749902/, 
-0.9973574665 ± 0.07265042297/. 

These agree numerically f2| with the eigenvalues of the monodromy matrix for y{s; 1). 

To get a better estimate of the value of m between 0.54 and 0.55 at which the orbit 
'Y{s;m) loses spectral stability as m decreases, we numerically computed Y{k/4-) for the 
values of m = 0.531,0.532, . . . ,0.538,0.539, and then computed the values of A-i and A,2. 
These show for m — 0.531 through m = 0.538 that 'Y{s;m) is linearly unstable because Xi 
and A,2 form a complex conjugate pair with nonzero imaginary part. For m — 0.539, we 
have that y{s;m) is linearly stable because 

A-i =0.1425261155, A.2 = 0.08595095311, 

which are real, distinct, have absolute value smaller than 1, and none are equal to or 
± 1 / ^/2. These eigenvalues of K imply that the characteristic multipliers of y{s; 0.539) are 
1 with algebraic multiplicity 4, and the two complex conjugate pairs 

0.8407916212 ±0.5413588917/, 
0.9413360780 ± 0.3374705738/, 

with each one of these having modulus 1. Thus the value of m in [0.53,0.54] at which 
'Y{s;m) is at least spectrally stable, lies in the interval (0.538,0.539). 
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